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Abstract The structures and distributions of light nuclei are investigated 
within a microscopic correlation model. Two particle correlations are re- 
sponsible for the scattering of model particles either to low momentum- or 
to high momentum-states. The low momentum states form the model space 
while the high momentum states are used to calculate the G-matrix. The 
three and higher order particle correlations do not play a role in the latter 
calculation especially if the correlations induced by the scattering operator 
are of sufficient short range. They modify however, via the long tail of the 
nuclear potential, the Slater determinant of the (A) particles by generating 
excited Slater's determinants. 
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1 Introduction 

Correlation effects in nuclei have first been introduced in nuclei by Villars pQ , 
who proposed the unitary-model operator (UMO) to construct effective op- 
erators. The method was implemented by Shakin [2j[3] for the calculation of 
the G-matrix from hard-core interactions. 

The UMO is based on the separation of the two body potential in a 
short and a long components. Within this separation the effective n-body 
Hamiltonian contains only the long component. The short-range component 
is considered up to the two body correlation and produces no energy shift in 
the pair state. 

Non perturbative approximations of the UMO have been recently applied 
to even nuclei in Ref. [4,5] which is treated here in more detail. The basics 
formulas of the Boson Dynamic Correlation Model (BDCM) presented in the 
above quoted paper have been obtained by solving the n-body problem in 
terms of the long range component of the two-body force. This component 
has the effect of generating a new correlated model space (effective space) 
which departs from the originally adopted one (shell model). The amplitudes 
of the model wave functions are calculated in terms of non linear equation 
of motions (EoM). 

By linearizing the systems of commutator equations, which characterize 
the EoM, we derive the eigenvalue equations for our model space. Within this 
correlated formalism we generate a model that includes not only the ladder 
diagrams of Ref. [6] but also the folded diagrams of Kuo [7] . 

The n-body matrix elements which define the eigenvalue equations are 
calculated exactly via the Cluster Factorization Theory (CFT) [5]. 

In this paper the BDCM model is applied to calculate the influence of 
the correlations on the charge distributions of the lithium isotopes and on 
the charge distribution of 6 He. 

The value obtained for the charge radius of the correlated 6 He is slightly 
bigger than the radius calculated in other theories [9-13] and that derived 
within the isotopic-shift IS theory [14]. A charge radius which agrees with 
the radii calculated in the Refs. [9-13] and those calculated in the cluster 
models of Refs. [15-17] is on the other hand obtained by considering only 
two protons in the si. This non correlated radius agrees also with the radius 
derived at Argonne within the IS theory [14] . Correlations have therefore the 
property to increase the charge radius of 6 He as observed for the isotopes of 
Lithium. 

The calculations performed in Ref. [18] for the charge radii of the lithium 
isotopes, although in good agreement with those measured at GSI-TRIUMF [19- 
21] and analyzed with the help of Ref. [22], are always slightly larger than 
those measured. For the stable isotope 6 Li the calculated radius agrees with 
the value obtained with the electron scattering experiments of Ref. [23] . How- 
ever, the charge radii calculated in the IS theory could also depend on the 
nuclear correlations. The calculations performed in Ref. [24] for the field shift 
(FS) of 7 Li show that the departure from a point nuclear approximation is a 
rather large effect. Additionally the higher order cross term contributions of 
Ref. [25] need to be considered. 
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A direct comparison between the calculated and the measured charge 
radii should therefore be performed after an accurate analysis of these two 
correcting factors. 

2 Theory of Two Correlated Particles 

In order to describe the structures and the distributions of nuclei we start 
from the following Hamiltonian: 

H = ^2(a\t\P) a) a a f3 + ($ a p\vi 2 \$i S ) a^a^agOy (1) 

a/3 a/?7^ 

were v%2 is the singular nucleon-nucleon two body potential. Since the ma- 
trix elements \uf3) are uncorrelated the matrix elements of v\2 are infinite. 
This problem can be avoided by taking matrix elements of the Hamiltonian 
between correlated states. In this paper the effect of correlation is introduced 
via the e lS method. In dealing with very short range correlations only the S2 
part of the correlation operator need to be considered. 

Following Ref. [2,3] we therefore calculate an "effective Hamiltonian" by 
using only the S2 correlation operator obtaining: 

Heff = e- iS *He iS > = £ a/J <a|i|/3)at 0/J + Ea^^l^l^^^ojo^ 

where v[ 2 is the long component of the two body interaction (note that the 
v[ 2 is hi the following equations simply denoted as v). The is the two 
particle correlated wave function: 

& a0 = e lS2 $ al3 (3) 

In dealing with complex nuclei however the (Si, i = 3 • • • n) correlations 
should also be considered. 

The evaluation of these diagrams is, due to the exponentially increasing 
number of terms, difficult in a perturbation theory. 

We note however that one way to overcome this problem is to work with 
e %(Si+S 2 +S 3 -\ hS n ) p era t or on the Slater's determinant by keeping the n- 
body Hamiltonian unvaried. 

After having performed the diagonalization of the n-body Hamilton's op- 
erator we can calculate the form of the effective Hamiltonian which, by now, 
includes correlation operators of complex order. 

Using Eq. |T]), we can compute the commutator of the Hamiltonian with 
the operator iflj^^Y that creates a valence particle pair. By performing 
this calculation we shall retain the linear and the non-linear terms which are 
formed by coupling the shell model states with the particle-hole excitations 
of the core. 

At this point in order to obtain a complete system of equations, we also 
have to calculate the commutator of the Hamiltonian |T]) with the non linear 
terms. With this calculation we introduce in the model space states which are 
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formed by coupling the valence states with the two particle-hole excitations 
of the core. 

The successive model equations are then formed by calculating the com- 
mutator with the operator obtained in the previous step. 

The set of commutator equations above is suitable to be solved by means 
of a perturbation expansion. The perturbative solution of the system of com- 
mutator equations is however not easily obtainable due to the high number 
of diagrams one needs to calculate. 

A non perturbative solution of the system of commutator equation can 
be obtained within the linearization method, which consists by applying the 
Wick's theorem to the ((n+2)p-2h) terms and by neglecting the normal order 
terms. 

This approximation is motivated by the consideration that the low lying 
spectra of nuclei the ((n+2)p-2h) terms are lying at much higher energy than 
that of the ((n+l)p-lh) states. 

The linearized system of the commutator equations is then solved exactly 
in terms of the CFT which calculates the n-body matrix elements in an 
expedite and exact way. 

In the following we give the basic formula of the method for a (nuclear) 
system with two valence particles. In second quantization, the two particle 
states are defined by: 

— = = [aia} 2 ] J M \0), (4) 

where, for the sake of simplicity, we have omitted the isospin quantum num- 
bers and where 

ai <-» hh (5) 

has been introduced to ensure a compact index notation of the angular mo- 
menta of the two particles. In this notation, the operator product a] i a] 2 just 
creates two coupled particles of single particle ji and j2 coupled to the final 
J quantum number. 

To derive the effect of the correlation on the two valence particles, we 
have, at this stage, to evaluate the next commutator 

[H,A\( ai J)]\0) 

= [£>aia Q + \ £ (aPWr^alalasa^ia^ym (6) 
which, after some operator algebra becomes 



[H,A\( ai J)]\0) 

= J2f2(2p\2 P ')A\((3 1 J)]\0) + Yl fl(2p|3pl/04(/MJ£J)]|0). (?) 

01 02 J[J' 2 
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In Eq. the A\(f3iJ) operators are those of Eq. g]) and the A\((3 2 J[ J 2 J) 
are defined below: 

$ 3p lk — 'I'yyyy = 4(&W)|0) = {{a), a),)^), a^Y \Q) . (8) 

J 1 J 2 J 3 J 4 J 1 J 2 J 3 

In Eq. ((U we have used the additional convention: 

P2 —> j'.h.iU\ (9) 

and we have associated: 

J{ to the coupling of j[ f 2 ^ 
J 2 to the coupling of j' 3 j' 4 ^ ' 

Having extended the commutator as in Eq. ([7]), we have also to calculate the 
commutator equation for the A\{a 2 J\ J 2 J) operators as given below: 

[H,Al{a 2 .hJ 2 J)]\Q) _^ 

M{Zplh\?>p'lh')A(l3 2 J' 1 J' 2 J)\Q))+ ^plh\ip2h)Al(/3 3 JiJ 2 JsJ)\0), 

(11) 



where we have introduced the (4p-2h) wave functions defined below: 

$A V 2h — > $ J } J '? J k J 'f J , = A\(f3 3 JUUU)\0) , x 

W n iK 1 2 d Jl ' 12) 

= (((a],at,)^(a],a i4 )^)^(4^) 4 ) J |0}, 
and where we have consistently extended the definition given in (|5ll0p : 

h ■ Aj'lHU'jl, (13) 

with: 

J[ associated to the coupling of j[ j' 2 

J' 2 associated to the coupling of j' 3 j' 4 (14) 

J3 associated to the coupling of j' 5 j' 6 

In the definition of A\(f} 3 J[ J 2 J 3 J) the coupling of J[ to J 2 to J' 12 has been 
discarded from the notation. In Eqs. (|7I11[) the J?'s are the matrix elements of 
the Hamilton's operator in the model wave functions. The next step would 
then be the computation of the commutator of the Hamiltonian with the 
A\{f3 3 J[ J 2 J3 J) operators. Here we linearize these contributions by consider- 
ing that in the study of the low energy spectrum and in the calculation of 
ground-state correlated distributions the A\(f3 3 J[ J 2 J 3 J) terms have a small 
contribution. The linearization is performed by applying to the (4p2h) terms: 

Y, (a0\v(r)\j5)aiala 5 a^At(p 3 J[4j' 3 J) (15) 
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the Wick's theorem and to discard the normal order terms. Within this lin- 
earization approximation we generate non perturbative solutions of the EoM 
from the commutator equations of Eq. (|7I11[) . i.e.: the eigenvalue equations 
for the mixed mode system: 

[H,A\{ ai J)]\0) = ^(2p|2p')4(M|0) 

131 (16) 

+ J2 n(Mzp'ih')4(M4J)\o)i 



and 



[H, A\(a 2 J 1 J 2 J)\\Q) = ^tf(3pl/qV)4(/? 1 J)|0) 

01 (17) 
+ ^{^h\Zp'lh')Al{p 2 J' x J' 2 )\Q). y ' 

foj[j' 2 

Within the application of the GLA approximation we convert Eqs. (ffl fTTj) in 
an eigenvalue equation for the configuration mixing wave functions (CMWFs) 
of the model. In fact, the linearization provides the additional matrix ele- 
ments necessary to write the following identity: 

Q^plhWlh') = hhhU '■:r:. l \j',i!,,'. . (18) 

and to introduce the off-diagonal matrix elements which couple the (2p) to 
the (3plh) subspaces. Now, by writing Eqs. (|16I17[) in the following matrix 
form: 

[H,A\( ai J)}\0) 

ff,4(a 2 JiJ a J)]|0>, . 
E 2p + Q(2p\2p') n(2p\3p'lh') \f A{(fi!J)\0) \ 1 ' 
Q(3plh\2p') E 3plh + f2(3plh\3p'lh') J \A\{p 2 J x J 2 J)|0) 



and by multiplying to the left with: 

(0|Ai(aiJ) 
(0\A 2 (a 2 JiJ 2 J) 

we generate the eigenvalue equation for the dressed particles: 



(20) 



(E 2p + {A l {a l J)\v{r)\A\{l3 l J)) {A x {a\ J) \v(r)\ A\ (f3 2 J[J 2 J)) 

R t1, v \ {A 2 (a 2 J 1 J 2 J)\v(r)\A{(p 1 J)} E 3plh + {A 2 {a 2 J x J 2 J)\v{r)\A\{(3 2 J[J' 2 J)) 

Xi{P\J) \ _ TP ( Xi{&iJ) 



X2W2W) )~ E [x2(a 2 J 1 J 2 J) ) |0) ' 

(21) 

In Eq.® E 2p = ef/ + ef/ and E 3plh = ef t F + ef/ + eg F - ef/ are the 
Hartre-Fock energies while the x' s are the projections of the model states: 

\H P ) = Xi(aiJ)^i(aiJ)|0) + X2 («2 Ji J2 J)A 2 (a 2 J t J 2 J)|0) (22) 
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to the basic vectors 2p, 3plh. To conclude, although the (4p-2h) CMWFs are 
not active part of the model space, they are important for structure calcula- 
tions. One may therefore associate the GLA approximation to a parameter 
which describes the degree of complexity of the model CMWFs. Within the 
first order linearization we obtain the EoM for the shell model while within 
the second and third order linearization approximations we derive the EoM 
of valence particles coexisting with the complex particle-hole structure of the 
excited states. 

In this paper we solve Eq. pTjl self-consistently. The solutions for the first 
iteration step are obtained by diagonalizing the eigenvalue equation (|2ip . The 
first step of the iterative method generates the dynamic amplitudes for the 
two dressed particles, i.e. two particles coexisting with the 3plh structures. 
With the calculated eigenvectors we recompute then the matrix elements 
(jlh\v(r)\jij2j334) and (jij2j3j4\v(r)\j' 1 j 2 j' 3 f 4 ) and we diagonalize again the 
eigenvalue equation. The iterations are repeated until the stabilization of the 
energies has been reached. Before performing the diagonalization of relative 
Hamilton's operator in the CMWFs base we have to eliminate the spurious 
center of mass components. In the BDCM this is performed, following the 
calculations of Refs. [26, 27], by calculating the percent weights of spurious 
states in the model wave functions. These can be obtained by evaluating the 
energy of the center of mass according to the following equation: 



(23) 



E R = J dR^ dressed { Jl3] J){R 2 )^ dressed { h jjJ) 



In Eq. (|23p the calculation of the integrals can be performed by using the 
CFT expansion for the (3plh) states and by considering that for two particle 
states we have: 



(j l j j J\(r i -r j )\j i j J J) 
- ±n\lj.] ( i\ 1 h 



where: 



J = (2.7 + 1). (25) 



By diagonalizing the above operator in the model space we obtain the energy 
of the center of mass. The overlap with the model space give the degree of 
"spuriosity" of the different components. The model space which characterize 
the BDCM is formed by adding even particle to a closed-shell nucleus. The 
closed shell configuration can be described by a single Slater determinant 
and one can use the Hartree-Fock's theory to obtain the binding energy and 
the single-particle energies. Alternatively one can remark that for a closed 
shell nucleus (Z,N) the single particle energies for the states above the Fermi 
surface are related to the binding energies differences: 

e> =BE(Z,N) -BE*(Z + l,N), (26) 
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Table 1 Single-particle scheme and single particle energies (MeV) used to form 
the model CMWFs for the A=6 isotopes 



hole 
energy 


lSl/2 

-20.58 
















hole/particle 
energy 


1^3/2 

1.43 
















particle 




1^5/2 


2si/ 2 


1^3/2 


l/V/2 


2P3/2 


1/5/2 


2P1/2 


energy 


1.73 


17.21 


22.23 


23.69 


25.23 


27.18 


28.33 


29.67 



and 

e > = BE{Z,N) -BE*{Z,N + 1). (27) 

The single particle energies for the states below the Fermi surface are given 
by: 

e< = BE*(Z-l,N)-BE(Z,N), (28) 

and 

e< =BE*(Z,N- 1) -BE(Z,N). (29) 

The BE are ground states binding energies which are taken as positive values, 
and e will be negative for bound states. (BE* = BE~E X ) is the ground state 
binding energy minus the excitation energy of the excited states associated 
with the single particle states. Within this method, which recently has been 
reintroduced by B.A. Brown [2S], we derive the single particle energies from 
the known spectra of neighbor nuclei (see Table (1)). 

The generalization of the previously defined formalism needed to calculate 
the charge radii of the A=7 (one particle DCM), A=8 (four particles BDCM) 
and A=ll (three particles DCM) is not given explicitly in this paper, but 
will be presented shortly. 

2.1 Results 

In order to perform structure calculations, we have to define a single particle 
base with the relative single-particle energies and to choose the nuclear two- 
body interactions. The single-particle energies of these levels are taken from 
the known experimental level spectra of the neighboring nuclei and given in 
Table (1). For the experimentally unknown single particle energies of the fp 
shells we use the corresponding energies for the mass A=9 nuclei scaled ac- 
cordingly the different binding energies. In this paper we perform as in [12,13] 
calculations by assuming all levels as bound for the particle-particle interac- 
tion, we use the G-matrix obtained from Yale potential [29) . These matrix 
elements are evaluated by applying the e correlation operator, truncated 
at the second order term of the expansion, to the harmonic oscillator base 
with size parameter b=1.76 fm. As also explained in Ref. [4,5] the potential 
used by the BDCM is separated in low and high momentum components. 
Therefore, the effective model matrix elements calculated within the present 
separation method and those calculated by Kuo [30-33] are pretty similar. 
The separation method generates matrix elements, which are almost inde- 
pendent from the radial shape of the different potentials generally used in 
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Table 2 Calculated charge radii for 6 He in fm compared with the results obtained 
in other theoretical models and with the radius derived within the IS theory. 



charge radius of He 



Model 



1.944 [9 no-core shell model 

2.09 [H quantum Monte Carlo technique 

2.25 this work BDCM 

2.39 this work BDCM without the folded diagrams 

2.06 this work two correlated Is i -protons 

1.99 ED Cluster' 

1.99 m Cluster 

1.99 Cluster 

2.054 ± .014[T4] Isotopic Shift (Exp.) 



structure calculations. 

The particle-hole matrix elements could be calculated from the particle- 
particle matrix elements via a re-coupling transformation. We prefer to use 
the phenomenological potential of Ref . [34] . The same size parameter as for 
the particle-particle matrix elements has been used. In Table (2) the cal- 
culated charge radii of 6 He are compared with the radii calculated by the 
other theoretical models and with the radius obtained by the IS theory. In 
the theoretical models quoted in this table the calculations of the charge 
distributions and of the charge radii are performed in terms of non- corre- 
lated operators. The correlations are included only in the derivation of the 
S2 effective Hamiltonian. 

For the stable 6 Li the calculated charge radius is equal to 2.55 fm, a value 
that reproduces well the charge radius of 2.55 fm obtained in Ref. [23] from 
the electron scattering experiments. The charge radii given in Table (3) for 
other lithium isotopes are however larger then those calculated in the other 
quoted theoretical models and then those obtained within the IS theory. 
Here also the main difference between the results obtained in the DCM and 
BDCM models and those of the other theoretical calculations has to be found 
in the treatment of the correlation operator. In [9,10] the charge radii are 
calculated in the "no core shell model" which is based on exact solutions of 
the two particles Schrodinger's equation by considering large computational 
spaces. The calculations do not include however the S3 correlations. [12-13] 
presents charge radii evaluated within an accurate Quantum Monte Carlo 
Method. The Hamiltonian used includes two- and three-bodies forces in a 
two-body correlated mechanics. The calculation method of [36-38] is based 
on a microscopic cluster method in which few particles arc interacting with 
the rest nucleus considered in its ground state. In our model the excitations 
of the core are associated to the £3 correlation operator which increases the 
charge radius of the Lithium isotopes. 



3 Conclusions and Outlook 



In this contribution we have investigated the effect of the microscopic correla- 
tion operators on the charge distributions of 6 He and of the Lithium isotopes. 
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Table 3 Calculated charge radii for the Lithium isotopes in fm compared with the 
results obtained in other theoretical models and with the radius derived within the 
IS theory (Exp.). 



Lithium 


Exp.(GSI) 19 


Exp.+Theo. [35J 


Theo. |9] 


Theo. 12 


Theo. 36 


DCM+BDCM 




rms 


rms 


rms 


rms 


rms 




b Li 


2.51 


2.47 


2.22 


2.54 




2.55 


7 Li 


2.39 


2.43 


2.13 


2.41 


2.43 


2.41 


8 Li 


2.30 


2.42 


2.13 


2.26 


2.34 


2.40 


9 Li 


2.22 


2.34 


2.16 


2.21 


2.27 


2.42 


n Li 


2.47 


3.01 






2.57 


2.67 



The microscopic correlation has been separated in short- and long-range cor- 
relations according the definition of Shakin [2,3]. The short-range correlation 
has been used to define the effective Hamiltonian of the model while the 
long-range is used to calculate the structures and the distributions of exotic 
nuclei. As given in the work of Shakin, only the two-body short-range cor- 
relation need to be considered in order to derive the effective Hamiltonian 
especially if the correlation is of very short range. For the long range corre- 
lation operator the three body component is important and should not be 
neglected. Within the three body correlation operator one introduces in the 
theory a three body interaction which compensates for the use of the genuine 
three body interaction of the no-core shell model. 

By using generalized linearization approximations and cluster factoriza- 
tion coefficients we can perform expedite and exact calculations. 

Within the calculated correlated distributions we obtain charge radii 
slightly larger than those calculated for non correlated distributions and de- 
rived by the IS experiments. 

The application of the DCM [24] to the two and three electron ener- 
gies and distributions of the Helium and Lithium atoms respectively could 
serve as future motivation for a reevaluation of the IS theory. From one side, 
the calculation of CM of the different isotopes, could help to obtain a non- 
perturbative formulation of the isotopic change of the electron transition 
energies. From the other side the field shift theory could include the correct 
isotopic variation. Since the derivation of the charge radii from the two pho- 
ton experiments is influenced by the precision of the theoretical calculations, 
this new proposed method could contribute to evaluate with better precision 
the charge radii of exotic nuclei. 
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